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Fermi Large Area Telescope (LAT) data analyses based on event reconstruction and classification 
are so far restricted to events of measured energy larger than 100 MeV. We present a new technique 
to recover the signal from Gamma-Ray Bursts' (GRB) prompt emission between ^30 MeV and 100 
MeV, which differs from the standard LAT analysis. Filling the gap between the energy ranges 
where the Gamma-ray Burst Monitor (GBM) and LAT operate is important to better constrain the 
high-energy spectra of GRBs. The LAT Low-Energy (LLE) technique is described, first performance 
studies are presented, as well as preliminary spectral re-analyses of two Fermi GRBs. 



I. INTRODUCTION 

Since the launch in June 2008, the Gamma-ray 
Burst Monitor (GBM) onboard Fermi detected over 
300 Gamma- Ray Bursts (GRB), 12 of which were de- 
tected and studied above 100 MeV with the Fermi 
Large Area Telescope (LAT). The GBM enables GRB 
spectroscopy from 8 keV to 40 MeV, and the standard 
LAT analysis is performed using events of energies 
above 100 MeV only. 

Here we propose a non-standard LAT analysis 
which allows to recover the prompt emission from 
GRB between -30 MeV and 100 MeV, joining the 
two instruments' energy ranges. This so-called LAT 
Low-Energy (LLE) technique also improves the pho- 
ton statistics above 100 MeV, and will allow us to 
better define the GRB prompt emission spectral char- 
acteristics. It is important to note that this analysis 
is based on non-public data and software, and is still 
under improvement and calibration at the time of this 
proceeding. 

Section [TT] describes the LLE technique objectives 
and principles. Peformances studies are presented in 
section lllll Preliminary re-analyses of the bright GRB 
080916C and GRB 090510 are presented in sectionQV] 



II. THE "LAT LOW-ENERGY" TECHNIQUE 

A. Motivation 

The LAT standard analysis is based on a sophisti- 
cated reconstruction and classification procedure [l|. 
The probability for each event to be a photon is esti- 
mated using the event topology and the reliability of 
its reconstruction. 
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FIG. 1: GRB 090510 spectral energy distribution 0. 
Top: the time-integrated spectrum clearly shows an ad- 
ditional power- law component. Bottom: time-resolved 
spectroscopy. In bin 'a' no signal was detected in the LAT. 
Bins 'b' and 'c' show an additional power-law component. 
No signal is detected in the GBM after 0.9 s (bin 'd'). The 
40 - 100 MeV gap is highlighted. 



The characteristcs of the particle backgrounds be- 
low 100 MeV is still under study and makes this clas- 
sification more complicated. Moreover, the LAT ef- 
fective area suffers from systematic effects which are 
not well determined yet in this energy range. Al- 
though the standard LAT analysis currently starts at 
100 MeV, public datasets and responses should allow 
for analyses starting at energies lower than 100 MeV 
for any kind of source in the future. 

On the one hand, faint or soft GRB emissions (e.g. 
bin 'a' of GRB 090510, see fig. [TJ are not detected so 
far above 100 MeV, while they can exhibit additional 
components which become dominant above few tens 
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FIG. 2: GRB 080916C prompt emission between 30 MeV 
and 100 MeV. Top: the lightcurve is used to define "ON" 
and "OFF" emissions. The "OFF" rate is extrapolated 
to estimate the background level in the "ON" interval. 
Bottom: signal and extrapolated background spectra in 
the "ON" interval. 



of MeV. Their signal above this energy can be rejected 
by the usual quality selections or be intrinsically too 
faint. An improved spectral analysis above 40 MeV 
would reveal any kind of feature, e.g. a faint addi- 
tional component, a cutoff, or a simple power-law be- 
haviour. The use of the present technique should im- 
prove the statistics above 100 MeV as well, although 
it is not yet clear whether it will substabtially improve 
the sensitivity of the spectral analysis. 

On the other hand, bright spectra (e.g. bin V to 
'd' of GRB 090510, see fig. [Q or GRB 080916C Q) 
are currently reconstructed above 100 MeV. Adding 
a LAT low-energy signal to the available standard 
dataset would help to better constrain the spectral 
parameters of such features (see section HV|) . 



B. Principle 

The LLE technique makes use of the analysis soft- 
ware RMfit, which is commonly used for all GBM 
spectral and temporal analyses. It consists of the 
forward- folding analysis of a background-subtracted 
binned event rate, using a Detector Response Matrix 
(DRM) for model- folding (fig. [3]). A temporal selec- 
tion defines the "ON" interval containing the emis- 
sion and the "OFF" interval(s) used to extrapolate 
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FIG. 3: Top: DRM produced for GRB 080916C LLE spec- 
tral analysis in the 10 MeV - 1 GeV energy range. Bot- 
tom: GRB 080916C counts spectrum between 30 MeV 
and 200 MeV (points) with folded spectrum (line) and 
residuals (bottom panel). 



the background rate in the "ON" interval, in each en- 
ergy bin (fig. [2]). 

In contrast with the standard analysis, the LLE 
event selection does not include strict quality crite- 
ria. We consider every event passing the onboard 
GAMMA filter [l| and for which at least one track 
could be found in the tracker. Many of these events 
do not have a reliable direction, therefore no "region of 
interest" can be defined around the source position in 
the sky. A binning is applied on the measured energy 
(see section UlI Bp . 

Once the "ON" and "OFF" intervals have been se- 
lected, the background rate time profile is fitted with 
a polynomial function in each energy bin, which yields 
an estimate of the background spectrum in the "ON" 
interval. The background spectrum can vary with 
time so the "OFF" intervals can not be chosen too 
far from the "ON" interval, nor too close nor short to 
allow for a good fit precision. The same precautions 
are usually taken for spectral analyses of GBM data. 

The DRM is built from a dedicated extensive simu- 
lation of a bright point source using the gleam soft- 
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FIG. 4: Fraction of events from a large set of simulated 
photons passing the standard Pass6 "transient class" qual- 
ity selections (squares) and LLE relaxed selections (cir- 
cles), both normalized to the onboard GAMMA filter. 



TABLE I: LLE energy measurement resolution for simu- 
lated photons of inclination angle < 40° (top) or 40° < 
< 70° (bottom). Energies are in MeV. True energy Emc 
lines of width 4 MeV were selected. The resolution is here 
defined as RMS/<E mes > and the error as (<E mes > - 
Emc)/Emc- 



Emc 




RMS 


resolution 


error 


30 


27 


10 


37% 


-10% 


50 


45 


16 


36% 


-10% 


100 


90 


27 


30% 


-10% 


500 


490 


70 


14% 


-2% 


30 


30 


14 


47% 


0% 


50 


44 


18 


41% 


-12% 


100 


85 


34 


40% 


-15% 


500 


470 


80 


17% 


-6% 



ware [4]. The observation conditions reproduce those 
of the GRB of interest. In particular, the inclination 
angle of the source and livetime should be the same 
since both affect the effective area. The LLE selec- 
tions are applied to the dataset produced, the remain- 
ing events are binned in true and measured energy. 
The resulting distribution is converted to area units 
(fig. [3]). The effect of the simulated spectral shape 
and of the true energy binning on the reconstruction 
quality are yet under investigation. 

III. PERFORMANCES STUDIES 
A. Cuts efficiencies 

The LLE selection consists of the onboard GAMMA 
filter and the identification of at least one track in 
the tracker. This latter criterion ensures a better cor- 
relation between the true and measured energy (see 
section IlII Bp . 

These relaxed selections greatly improve the photon 
statistics at low energies with respect to the standard 
Pass6 "transient class" quality selections, and improve 
the statistics at high energies as well, as shown in fig. HI 
: the effective area is increased by a factor of ~50 at 
30 MeV, 4.5 at 100 MeV, and still by more than a 
factor 2 above 1 GeV. 



B. Energy resolution 

The energy estimate used for this analysis is dif- 
ferent from the standard one. It is based on energy 
depositions in both the tracker and the calorimeter. 
For low-energy events (below 100 MeV) which barely 
reach the calorimeter, the tracker mainly contributes 
to the energy estimate. Hence the requirement of at 
least one reconstructed track. The choice of this en- 
ergy estimate may not be final, though it nicely cor- 



relates to the true energy for simulated datasets (see 

fig. ED. 

A study of the energy resolution from this estimate 
has been performed, using large samples of simulated 
photons (see table I]). No very large bias is found, 
and the energy resolution is reasonable even at low 
energies (50% above <30 MeV, <40% above 50 MeV). 

C. Reconstruction capabilities 

To study the reconstruction capabilities of the LLE 
technique, a bright GRB was simulated as a point 
source with a simple power-law spectrum : N(E) — 
No(E/Eo)~P . No background was added to the data 
and this spectrum was analyzed using the standard 
tool XSPEC (HEASARC). The reconstruction was 
performed using only the standard Pass6 transient 
data (above 100 MeV), or using both the LLE (30 
MeV to 100 MeV) and Pass6 transient data. Both 
reconstructions yield similar results (table HI]) . with a 
smaller uncertainty on the index j3 when using LLE 
data. A systematic error yet appears, which will be 
investigated after the technique has been fully cali- 
brated. 



TABLE II: Spectral parameters of the simulated GRB 
spectrum, and reconstructed spectra using standard Pass6 
transient events only, or LLE events as well. Only statisti- 
cal errors are shown. Though the use of LLE data reduces 
the uncertainties, a systematic error appears. 



Spectrum 


iVo 

(10" 9 ph.cm^.s- 1 ) 


P 


input 


1.19 


2.1 


fit without LLE 
fit with LLE 


1.36 ± 0.13 
1.60 ± 0.12 


2.21 ± 0.06 
2.02 ± 0.02 
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FIG. 5: Counts spectra of a simulated bright GRB with 
folded model and ratio of the signal to the folded model. 
Top: only standard Pass6 transient events above 100 MeV 
are used, events converting at the front and back of the 
tracker are separated 5]. Bottom: LLE events between 
30 MeV and 100 MeV are added to the previous dataset. 
The ratio panel shows a good agreement between both 
sets, and the fit is dominated by the LLE data. 



IV. PRELIMINARY RE-ANALYSES 
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FIG. 6: Top: GRB 080916C counts spectrum with folded 
model and residuals. GBM and LAT standard transient 
data were fitted. LLE dataset was superimposed to the 
spectrum. Bottom: GRB 090510 time-integrated counts 
spectrum with folded model and residuals. GBM data, 
LAT Pass6 transient data (above 100 MeV) and LLE data 
(30 MeV - 100 MeV) were fitted together by a Band func- 
tion with an additional power- law component. 



New analyses of two bright bursts were performed 
and yield encouraging results. 

GRB 080916C time-integrated spectrum was fitted 
using only GBM and standard transient events (above 
100 MeV) (fig. Ej). The LLE dataset was added but 
not used for the fit. The good residuals in the LLE 
energy range (30 MeV - 100 MeV) show the good 
agreement between this technique and the standard 
procedure. 

GRB 090510 time-integrated spectrum was fitted 
using all datasets : GBM, LAT Pass6 transient, LLE. 
A Band function with an additional power-law com- 
ponent yields the best fit, like in the standard analy- 
sis [2| . The high-energy component is more significant 
when using LLE data (N a = 8.9 instead of 5.6 without 
LLE). The spectral evolution observed in this bright 
burst's prompt emission will certainly better show up 
if the LLE data are used in the time-resolved spec- 
troscopy. 



V. CONCLUSION 

The LLE technique presented here can be used in 
principle for any kind of flaring source : GRB prompt 
emission, AXP, pulsars (see [6[ for details), etc. 

In particular it appears to be very promising for 
GRB prompt emission spectral analyses, revealing or 
better defining GRB spectral features above 30 MeV. 

The validation study of these analyses is yet still on- 
going (acceptance and energy calibration, systematics 
effects). The reconstruction of a simulated spectrum 
has shown a bias, such effects have to be character- 
ized and understood. The performances that are yet 
measured from sets of simulated photons have to be 
confronted to data, e.g. using Vela on-pulse emission 
as a pure sample of real photons. 
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